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In cluster physics a single particle potential to determine the microscopic part of the total energy 
of a collective configuration is necessary to calculate the shell- and pairing effects. In this paper 
we investigate the properties of the Riesz fractional integrals and compare their properties with 
the standard Coulomb and Yukawa potentials commonly used. It is demonstrated, that Riesz 
potentials may serve as a promising extension of standard potentials and may be reckoned as a 
smooth transition from Coulomb to Yukawa like potentials, depending of the fractional parameter a. 
For the macroscopic part of the total energy the Riesz potentials treat the Coulomb-, symmetry- and 
pairing-contributions from a generalized point of view, since they turn out to be similar realizations 
of the same fractional integral at distinct a values. 
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I. INTRODUCTION 

Convolution integrals of the type 

F{x)^ I dUix-O^iO 
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m 
m 
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play a significant role in the areas of signal- and image 
processing or in the solution of differential equations. 

In classical physics the first contact with the 3D- 
generalization of a convolution integral occurres within 
the framework of gravitation and electromagnetic theory 
respectively in terms of a volume integral to determine 
the potential of a given charge density distribution p: 
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where the weight w 
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is interpreted as the gravitational or electromagnetic field 
of a point charge [221124]. 

In nuclear physics collective phenomena like fission or 
cluster-radioactivity, where many nucleons are involved, 
are successfully described introducing the concept of a 
collective single-particle potential, based on folded poten- 
tials of e.g. Woods-Saxon type. Several weight functions 
have been investigated in the past. 

In this paper we will demonstrate, that the fractional 
Riesz-potential [40l [HI |46] which extends the weight 
function introducing the fractional parameter a 
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serves as a serious alternative for commonly used Nilsson- 
[38l|39], Woods-Saxon [8] and folded Yukawa potentials 
[TJ 133) , modeling the single particle potential widely ap- 
plied in nuclear physics as well as in electronic cluster 
physics. 

Hence we give a direct physical interpretation of a 
multi-dimensional fractional integral within the frame- 
work of fragmentation theory |31j . which is the funda- 
mental tool to describe the dynamic development of clus- 
ters in nuclear and atomic physics. 



II. FOLDED POTENTIALS IN 
FRAGMENTATION THEORY 

The use of collective models for a description of collec- 
tive aspects of nuclear motion has proven considerably 
successful during the past decades. 

Calculating life-times of heavy nuclei [HI [371 155] . 
fission yields [ 28] . giving insight into phenomena like 
cluster-radioactivity [44], bimodal fission [TQ or mod- 
eling the ground state properties of triaxial nuclei [35] 
- remarkable results have been achieved by introducing 
an appropriate set of collective coordinates, like length, 
deformation, neck or mass-asymmetry [30] for a given 
nuclear shape and investigating its dynamic properties. 

As an example, in figure [Tj the parametrization of the 
3-spheres-model is sketched. It determines the geometry 
of a given cluster shape by 2 intersecting spheres, which 
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are smoothly connected via a third sphere, which models 
a neck depending on the size of the radius r^. 

The corresponding set of collective coordinates {q^,i = 
1, ...,4} is given by [7]: 



the two center distance: 

Az = 



Z2 - Zl 



the mass asymmetry: 



VA = 



A1-A2 
A1+A2 



the charge asymmetry: 



Zl — Z2 



Zl 



the neck: 
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where Ai, Z\ and A2, Z2 are the number of nucleons and 
protons in the two daughter nuclei. 

This choice of collective coordinates allows to describe 
a wide range of nuclear shapes involved in collective phe- 
nomena from a generalized point of view |17j . e.g. a si- 
multaneous description is made possible of general fission 
properties and the cluster-radioactive decay of radium 



223Ra ^ 209pb 
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which was predicted by Sandulescu, Poenaru and Grei- 
ner in 1980 and later experimentally verified by Rose and 
Jones in 1984 gSllSnj. 

In order to describe the properties and dynamics of 
such a process, we start with the classical Hamiltonian 
function 



= T + Vb 



(11) 




FIG. 1. The parametrization of the 3-spheres-model. 



introducing a collective potential Vq, depending on the 
collective coordinates, 



(12) 



with a macroscopic contribution i?macro based on e.g. the 
liquid drop model and a microscopic contribution E'mic, 
which mainly contains the shell and pairing energy based 
on a single particle potential Vs.^. and the classical kinetic 
energy T 



(13) 



with collective mass parameters Bij. 

There are several common methods to generate the 
collective mass parameters B^^ e.g. the cranking model 
[5U] or irrotational flow models are used |23) . 

Quantization of the classical Hamiltonian [42] results in 
the collective Schrodinger equation 



^o*(g\<)= - 



1 
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with B ~ det Bij is the determinant of the mass ten- 
sor. This is the central starting point for a discussion of 
nuclear collective phenomena. 

For a specific realization of the single particle poten- 
tial Vs.p., for protons and neutrons respectively a Woods- 
Saxon type potential may be used. The advantages 
of such a potential are a finite potential depth and a 
given surface thickness. Furthermore arbitrary geometric 
shapes may be treated similarly by a folding procedure, 
which yields smooth potential values for such shapes. 

For the three-spheres model, in order to define a cor- 
responding potential, a Yukawa-function is folded with a 
given volume, which is uniquely determined within the 
model: 



VV(r) 



Vo 



dV ^P (15) 



r — r \ a 



with the parameters potential depth Vq and surface thick- 
ness a. 

For protons, in addition the Coulomb-potential has to 
be considered, which is given for a constant density po 



1 



a Jv r\/a 
where the charge density is given by 

Ze 



(16) 
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Both potentials may be written as general convolutions 
in of type: 



V'typo(r) = C'typc / d^r p{r )u;typo(|r - 



(18) 
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with the weights 



wvid) = 
wc{d) = 



exp" 



-dj a 
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where d = |a; — ^| is a measure of distance on and a 
density, which is constant inside the nucleus 



p{r) 



Pq r inside the nucleus 
r outside the nucleus 



(21) 



Therefore the single particle potential Vs. p. is given by 

Vs.p.=VY + {]^+t-^Vc + na{VVYXp) (22) 

with is the eigenvalue of the isospin operator with + 1 
for protons and — ^ for neutrons, which guarantees that 
the Coulomb potential Vc only acts on protons. The last 
term is the spin-orbit term with the Pauli-matrices a, p\s 
the momentum operator and the strength is parametrized 
with K. This term is necessary to split up the degener- 
acy of energy levels with different angular momentum 
and to generate the experimentally observed magic shell 
closures. 

In the original Nilsson oscillator potential an addi- 
tional I?' term was necessary to lower the higher angu- 
lar momentum levels in agreement with experiment. For 
Woods-Saxon type potentials such a term is not neces- 
sary. Whether Riesz potentials are a realistic alternative, 
will be investigated in the next section. 

The solutions of the single particle Schrodinger equa- 
tion with the potential Vs.p. yield the single particle en- 
ergy levels, which are used to calculate the microscopic 
part of the total potential energy and contains two major 
parts, the shell and pairing corrections. 



En,ic{q') = £^shcu(g') + £^pair(g*) 



(23) 



III. THE RIESZ POTENTIAL AS SMOOTH 
TRANSITION BETWEEN COULOMB AND 
FOLDED YUKAWA POTENTIAL 



If we reinterpret the Riesz potential 

Vrz (r) = Crz I d\'p{f )wRz{\r-r\) 



with the weight 

vjRz{d) 



(d/a) 



< a < 3 



(24) 



(25) 



as the 3D- version of the general Riesz integral ([IS]) ap- 
plied to a scalar function p{f), we may treat and interpret 
the Coulomb {a = 1), Riesz- and Yukawa potentials sim- 
ilarly from a generalized point of view. 
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FIG. 2. For a spherical assumed shape (here 

^^^U) the potential for different weight functions is drawn. 
From top to bottom: Coulomb (a = 1.00), Riesz (a = 
1.25,1.75,2.25,2.50) and Yukawa (with a = 0.9[fm] weight. 
In order to compare the plots with cylinder symmetric shapes, 
potential is drawn in cylinder coordinates {z, p) for a sequence 
of p = 0.00, 0.25, 0.50, 0.75, 1.00 x Rq. J?o(^^^U) = 8.26[/m] 



In the following we will investigate the behavior of the 
Riesz potential with varying a, and compare its proper- 
ties with the cases of Coulomb and Yukawa weight func- 
tions. In a way, the parameter a in the Riesz potential 
may be interpreted as a global screening of the Coulomb 
weight, such that the effect of the Yukawa exponential is 
partly modeled. 



wc{d) 
WRz{d) 
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(d/a)"-i d/a 

W(d)=exp-'^/''-^ 
d/a 



(26) 
(27) 
(28) 



Therefore the Riesz potential could be an interesting al- 
ternative to the Yukawa potential in the case a ^ 1. In 
a way, we expect the screening properties of the Riesz 
potential for increasing a to result in an interpolation 
between Coulomb and Yukawa limit. 

Hence the fragmentation potentials used in a dynamic 
description of fission or cluster emission processes are an 
ideal framework to discuss and understand the properties 
of the Riesz integral. 



The integral (18) with the weights (26l-(28) may be 
evaluated analytically for a spherical nucleus with radius 



p{r) = poH{Ro - r) 
with the Heaviside step function H. 



(29) 
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FIG. 3. For the configuration ^^^U Pb Ne tiie 

Coulomb, Riesz (a = 2.50) and Yukawa (a = 0.9[fm]) is plot- 
ted for = 0.75iio (left column) and = 1.375iio (right 
column). 7?o(^^^U) = 8.26[/m]. 



In this case we have: 



v; 



iphcrc / 



typo 



Ro 



{r) = CtypcPo r^dr / sin{e )de x 



2-rr 



with 



\r — r \ = \ r 



+ r'2 — 2rr' cos{6') 



With the substitution u 



\J r'^ + r'"^ — 2rr' coa{9') 



(30) 



(31) 



(32) 



we end up with a double integral for spherical shapes: 



sphere 



type 



(r) = 27rCtypo Po dr r /rx 



Ro 



\/ (r+r')2 



^ (r-r')2 

27rpo 



du uw typciu) 
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L^typo / , ' ' 

' dr r X 



r-fr 



|r— r I 



duuwtypciu) 



(34) 



This integral is valid for any analytic weight w{u) and 
may be easily solved for the Coulomb, Riesz and Yukawa 



weight functions. We obtain: 



y^ph"-(^) ^ aC'c 



Ze/3 _ ' 
-Ro 2 2Rl' 



r<Ro 



(35) 



Vlt"%r) = Ana'^Cnz-, rr^ r.- x 

(a — 2) (a — 6)[a — 4) r 

f(r + i?o)3-"(r-(3-a)i?o)+ 
(i?o-r)3-"(r + (3-a)i?o) r < i?o 
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r>Ro 
1 



(36) 



(r + i?o)3-"(r-(3-a)i?o)- 



r>Ro 

r(l_(l + ^)e-«o/asinM^) r<i?o 



^^phere^^^ ^ ATra^'Cy X 



(37) 



l^(f^cosh(^)-sinh(i^)) r>Ro 

In figure [2] a sequence of these potentials is plotted for 
a spherical nucleus, ranging from Coulomb (a = 1.00) 
and Riesz potential with increasing a up to the Yukawa 
potential with parameter settings according to [1]. 

For large a > 2.00 the Riesz potential as well as the 
Yukawa potential model a finite surface thickness. 

A remarkable difference between both potentials fol- 
lows for small z. In this area, the Yukawa potential 
models a more Woods-Saxon type potential, while the 
Riesz potential may be compared with an harmonic os- 
cillator potential. But this behavior is restricted only to 
the lowest energy levels; for realistic calculations the en- 
ergy region near the Fermi- level is much more relevant. 
In this region, both potential types show a similar be- 
havior for 2.00 <a< 2.50. 

Hence both potentials seem interesting candidates for 
generation of realistic singe particle energy levels. 



For cylinder symmetric configurations the integral ( 18 ) 



cannot be solved analytically. Instead we switch to cylin- 
der coordinates {p, z, cj)}. With the distance 



dcyi = \J p"^ + p"^ - "^.pp cos{(t>' ) + {z-z'Y (38) 
wc have to solve the integral 



Vtype(p, Z) = Ctype / dp p dz d4} p{p , Z )wtypc(dcyl) 

Jv 

(39) 

numerically. 

In figure [3] we have solved ( 39 ) and compare the three 
different weights for the strong asymmetric cluster decay 



232 



U 



,208 



Pb 



24 



Ne 



(40) 



The Riesz potential allows for a smooth transition be- 
tween the Coulomb case and the Yukawa limit by varying 
a. Hence we obtain a direct geometric interpretation of 
the fractional parameter a. 

Up to now, we discussed the properties of the single 
particle potential, which is the starting point for a calcu- 
lation of the microscopic part of the collective potential. 
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On a macroscopic level the self energy of a given config- 
uration contributes to the macroscopic part of the nuclear 
potential as the Coulomb- and surface or more sophisti- 
cated Yukawa energy term in a macroscopic energy for- 
mula, historically first used in Weizsackers famous liquid 
drop mass formula [60]: 

a,y^{N - Z)A-' + (41) 

as a function of the nucleon number A — toRq containing 
a volume, surface. Coulomb, symmetry and pairing term. 
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FIG. 4. For a spherical shape the self energy as a function 
of the sphere radius -Roffm] is plotted for the Coulomb, the 
Riesz (a = 1.49, 1.98, 2.47, 2.96 and the Yukawa (a = 0.9[/m] 
is plotted. To compare all different types, all energies are 
normalized to i5typo(iio ~ 1) ~ 1- Depending on Ro, the 
Yukawa self energy lies is covered by the Riesz self energy 
within a range of 2.0 < a2.5 values 



{Ze = 1) and with a thickness parameter a > 
3 a 



En — 



5Rn 



E 



9 X 22-"a" 
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RZ 



Ey = 



(3-a)(4-a)(6-a) R^ 
3a3 - 3aRl + 2i?g - 3ae-2«o/''(a + Rof 



(44) 

< a < 3 (45) 
(46) 



We obtain the important result that the Riesz self en- 
ergy behaves like 



Erz 



1 



which scales with the nucleon number ^ ~ i?n as 



E 



RZ 



(47) 



(48) 



and therefore allows to model the influence of a screened 
Coulomb like charge contribution to the total energy. 

In figure |4] we compare the i?o dependence for the 3 
different types of self energy. Depending on the size of 
the spherical nucleus e.g. Ro{^°^Pb) = 7.25 the behav- 
ior of the Yukawa self energy is covered by the Riesz self 
energy for a within the range 2 < a < 2.5 and there- 
fore the Riesz self energy covers the full range of relevant 
categories. 

In addition it should be mentioned, that for the case 
a = 3/2 the Riesz self energy behaves like A~^/^ which 
emulates the pairing term and for a = 3 the Riesz self en- 
ergy behaves like A~^ which is equivalent to the behavior 
of the proton-neutron symmetry term in the Weizsacker 
mass formula. 

As a consequence, the pairing-, symmetry- and 
Coulomb contributions to the total energy content of a 
nucleus may be treated from a generalized view as differ- 
ent realizations of the same Riesz potential. 



IV. CONCLUSION 



The self energy for a given charge type i^type is defined 
as the volume integral over the potential restricted on 
the volume of a given shape. 



;^Ctypo / d^rp{r) / d^r' p{r)wtype{\f~r\) 
^ Jv JR^ 

(42) 



Etypo{r) = ^Ctypo j j ^ d^rd^r p{r)p{r)wtyp^{\f-r\) 

(43) 



For the 3 different weights (26l-(28) we obtain for the 



simplest case of a sphere with rac 



ius i?o , a unit charge 



From all these presented results we may draw the 
conclusion, that the Riesz potential may be considered 
as a promising alternative approach to folded poten- 
tials, which are widely used to describe nuclear dynamics 
within the framework of a collective shell model. 

Of course, these potentials are only an alternative 
starting point to calculate fragmentation potentials based 
on a fractional integral definition, but it is indeed remark- 
able, that the Coulomb-, pairing- and symmetry part 
of the macroscopic energy contribution may be consid- 
ered as specific realizations of the fractional Riesz inte- 
gral with the fractional parameters a € {1,3/2,3} . 
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